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ABSTRACT 

The  equilibrium  structure  ot  homogeneous  turbulent  shear  flow  is  investigated  from  a 
theoretical  standpoint.  Existing  turbulence  models,  in  apparent  agreement  with  physical 
and  numerical  experiments,  predict  an  unbounded  exponential  time  growth  of  the  tur¬ 
bulent  kinetic  energy  and  dissipation  rate;  only  the  anisotropy  tensor  and  turbulent  time 
scale  reach  a  structural  equilibrium.  It  is  shown  that  if  vortex  stretching  is  accounted  for 
in  the  dissipation  rate  transport  equation,  then  there  can  exist  equilibrium  solutions,  with 
bounded  energy  states,  where  the  turbulence  production  is  balanced  by  its  dissipation. 
Illustrative  calculations  are  presented  for  a.  k  —  e  model  modified  to  account  for  vortex 
stretching.  The  calculations  indicate  an  initial  exponential  time  growth  of  the  turbulent 
kinetic  energy  and  dissipation  rate  for  elapsed  times  that  are  as  large  as  those  considered 
in  any  of  the  previously  conducted  physical  or  numerical  experiments  on  homogeneous 
shear  flow.  However,  vortex  stretching  eventually  takes  over  and  forces  a  production- 
equals-dissipation  equilibrium  with  bounded  energy  states.  The  validity  of  this  result 
is  further  supported  by  an  independent  theoretical  argument.  It  is  concluded  that  the 
generally  accepted  structural  equilibrium  for  homogeneous  shear  flow  with  unbounded 
component  energies  is  in  need  of  re-examination. 


’Research  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Contract 
No.  NASl-18605  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in  Science 
and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665-5225. 

i 


1.  INTRODUCTION 


Homogeneous  turbulent  shear  flow  has  been  the  subject  of  a  variety  of  experimental, 
computational  and  theoretical  studies  during  the  past  four  decades.  The  popularity  of 
this  flow  lies  in  the  fact  that  it  accounts  for  an  important  physical  effect  —  the  alteration 
of  the  turbulence  structure  by  shear  —  in  a  simplified  setting  unencumbered  by  such 
complications  as  rigid  boundaries  and  mean  turbulent  diffusion.  Von  Karman^  first 
proposed  the  problem  of  homogeneous  shear  flow  which  gave  rise  to  some  mathematical 
studies  during  the  1950’s  (c.f.  Townsend^  and  Hinze^  for  a  review).  It  is  a  difficult 
flow  to  simulate  experimentally  and  the  first  tangible  suggestion  for  its  generation  was 
put  forth  by  Corrsin^  The  first  successful  experimental  realization  of  homogeneous 
shear  flow  in  the  laboratory  was  achieved  by  Rose®  and  was  followed  by  a  series  of 
landmark  experiments  by  Corrsin  and  co- workers®"®.  More  recently,  Tavoularis  and 
Karnik®  and  Rohr  et  al.^°  performed  more  exhaustive  measurements  of  homogeneous 
shear  flow  shedding  new  light  on  its  basic  structure. 

With  the  dramatic  increase  in  computer  capacity  achieved  by  the  late  1970’s,  direct 
numerical  simulations  of  homogeneous  turbulent  flows  became  possible.  Rogallo“  was 
the  first  to  conduct  a  direct  simulation  of  homogeneous  shear  flow;  the  results  that  he 
obtained  were  well  within  the  range  of  the  experiments  of  Tavoularis  and  Corrsin®.  Sub¬ 
sequently,  Bardina  et  al.^^  performed  a  coarse-grid  large-eddy  simulation  of  homogeneous 
shear  flow  and  Rogers,  Moin  and  Reynolds^®  conducted  fine-grid  128  x  128  x  128  direct 
simulations  which  further  clarified  its  structure. 


The  early  experiments  of  Rose®  and  Champagne,  et  al.®,  which  were  conducted  for 
relatively  weak  shear  rates  and  small  elapsed  times,  seemed  to  indicate  that  the  Reynolds 
stresses  (and,  hence,  the  turbulent  kinetic  energy)  asymptoted  to  equilibrium  values.  This 
asymptotic  state  is  consistent  with  the  production-equals-dissipation  equilibrium  that 
had  been  hypothesized  by  Townsend®  several  years  earlier.  However,  the  integral  length 


scales  were  still  growing  monotonically  at  the  end  of  the  Rose®  and  Champagne  et  al.® 
experiments  suggesting  the  strong  possibility  that  an  asymptotic  state  had  not  yet  been 
reached.  Subsequent  physical  experiments®"^®  and  direct  numerical  simulations^^"^®,  con-  ® 
ducted  for  stronger  shear  rates  and  larger  elapsed  times,  confirmed  this.  These  physical 
and  numerical  experiments,  along  with  alternative  theoretical  analyses  (see  Rogallo^^  and  “ 
Tavoularis^'*),  have  led  to  the  following  widely  accepted  physical  picture  of  homogeneous 
shear  flow:  _ _ 
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(1)  the  turbulent  kinetic  energj’  k  and  dissipation  rate  e  grow  exponentially  in  time 
at  the  same  rate, 

(2)  the  anisotropy  tensor  6,^  and  the  dimensionless  turbulent  time  scale  Sk/e  reach 
equilibrium  values  that  are  relatively  independent  of  the  initial  conditions  and  the  shear 
rate  S. 

It  should  be  noted  that  Speziale  and  MacGiolla  Mhuiris^®  have  recently  shown  that 
virtually  all  of  the  commonly  used  two-equation  turbulence  models  and  second-order 
closures  are  consistent  with  this  hypothetical  picture  of  homogeneous  shear  flow. 

An  alternative  physical  picture  of  homogeneous  shear  flow  is  presented  in  this  paper 
which  is  consistent  with  the  previously  conducted  physical  and  numerical  experiments, 
yet  at  the  same  time  is  more  intuitively  satisfying  in  that  it  excludes  the  occurrence  of 
unbounded  energy  growth.  In  particular,  it  will  be  shown  that  when  vortex  stretching  is 
accounted  for  in  the  dissipation  rate  transport  equation,  a  production-equals-dissipation 
equilibrium  results  in  which  the  turbulent  kinetic  energy  and  dissipation  rate  eventually 
asymptote  to  bounded  values.  Illustrative  calculations  are  presented  for  a  A)  —  e  model 
suitably  modified  to  account  for  vortex  stretching.  Consistent  with  physical  and  numeri¬ 
cal  experiments,  these  calculations  indicate  an  exponential  time  growth  of  the  turbulent 
kinetic  energy  and  dissipation  rate  for  St  <  30.  This  is  the  largest  elapsed  time  considered 
in  any  of  these  previous  experimental  studies.  However,  for  St  >  30,  vortex  stretching 
eventually  takes  over  causing  the  system  to  saturate  and  attain  an  equilibrium  structure 
with  bounded  kinetic  energy  and  dissipation.  In  the  sections  to  follow,  a  detailed  case  is 
made  to  establish  that  this  alternative  equilibrium  structure  of  homogeneous  shear  flow 
is  a  serious  possibility  that  could  have  important  implications  for  turbulence  modeling 
at  high  Reynolds  numbers. 

2.  THEORETICAL  BACKGROUND 


We  will  consider  the  incompressible  and  isothermal  turbulent  flow  of  a  viscous  fluid. 
The  governing  field  equations  are  the  Navier-Stokes  and  continuity  equations  given  by 


dvt  dP 


5x, 


=  0 


(1) 

(2) 
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where  Vi  is  the  velocity  vector,  P  is  the  modified  pressure  and  i/  is  the  kinematic  viscosity 
of  the  fluid.  As  in  the  usual  treatments  of  turbulence,  the  velocity  and  pressure  will  be 
decomposed  into  ensemble  mean  and  fluctuating  parts,  respectively; 


u,  =  u,  +  u,,  P  =  P  +  p.  (3) 

For  homogeneous  shear  flow,  the  mean  velocity-gradient  tensor  takes  the  form 

/  0  5  0  \ 

^.  0  0  0.  (4) 
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In  direct  numerical  simulations  of  homogeneous  shear  flow,  an  initially  isotropic  turbu¬ 
lence  is  subjected  to  the  constant  shear  rate  S  and  its  time  evolution  is  then  computed.  In 
laboratory  experiments,  an  initially  decaying  isotropic  turbulence,  created  downstream 
of  a  grid,  is  subjected  to  a  uniform  shear  rate  as  it  evolves  spatially.  The  two  problems 
are  related,  in  an  approximate  sense,  by  the  Galilean  transformation 

X  =  Xo  +  Uct  (5) 

where  Uc  is  a  characteristic  mean  velocity  that  is  typically  taken  to  be  the  centerline 
mean  velocity  of  the  uniform  shear  [hence,  dimensionless  time  St  =  S{x  —  x<,)/I7c]. 
For  the  remainder  of  this  paper,  we  will  only  consider  the  temporally  evolving  version 
of  homogeneous  shear  flow,  since  it  is  the  only  version  of  this  problem  that  is  exactly 
homogeneous. 

In  any  homogeneous  turbulent  flow,  the  exact  transport  equations  for  the  turbulent 
kinetic  energy  k  =  and  dissipation  rate  e  =  take  the  form^® 

ic  =  V  -  e  (6) 


C  —  Vcg  +  Pey  ~  ^ 


where 


^  dv-  ^  ^  _ diTi 
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are,  respectively,  the  production  of  turbulent  kinetic  energy,  the  production  of  dissipation 
by  mean  strains,  the  production  of  dissipation  by  vortex  stretching,  and  the  destruction 
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of  dissipation.  In  (8)  -  (9),  =  u~Uj  is  the  Reynolds  stress  tensor  and  a;  =  V  x  u  is  the 

fluctuating  vorticity  vector.  At  this  stage,  we  introduce  the  anisotropy  tensor  bij  defined 
as 


_  (t.j  -  f  M.j) 
2k 


(10) 


which  will  be  useful  in  the  analysis  to  follow. 

If  we  non-dixnensionalize  the  turbulent  kinetic  energy,  dissipation  rate,  and  time  as 
follows: 


where  kg  and  Cq  denote  initial  values,  then  the  exact  transport  equation  (6)  for  k  can  be 
rewritten  in  the  dimensionless  form 


k*  =  2  -  l')  b^2k*  (11) 

since  "P  =  —T12S  in  homogeneous  shear  flow.  It  should  be  noted  that  the  ratio  of 
production  to  dissipation  "P/e  can  be  related  to  bi2  through  the  identity 


^  (12) 

e  e 

Consequently,  if  any  two  of  the  quantities  P/e,  612  and  Sk/e  are  known,  the  remaining 
one  can  be  computed  using  (12).  Furthermore,  since  for  homogeneous  shear  flow  P/e  >  0 
and  Sk/e  >  0,  it  follows  from  (12)  that  bu  <  0.  Physical  and  numerical  experiments 
have  tended  to  indicate  that  bi2  and  Sk/e  reach  equilibrium  values  that  are  relatively 
independent  of  the  initial  conditions.  More  specifically,  these  experiments  suggest  the 
results 

«6.0,  (612)00^-0.15,  (-)  «1.8  (13) 

where  (  .  )oo  is  the  equilibrium  value  obtained  in  the  limit  as  t  —*  00.  From  (11)  and  (12) 
it  follows  that  for  P  1 

A:*  ~  exp(At*)  (14) 

where  the  growth  rate  A  is  given  by 


A  =  2 


(612)00 


0.14. 


(15) 


Equations  (13)  and  (14)  then  imply  that  for  P  1 


e*  ~  exp(At*) 


(16) 
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since,  if  k  grows  exponentially  and  Skje  equilibrates,  then  e  must  grow  exponentially  at 
the  same  rate.  This  constitutes  an  alternate  derivation  of  the  Tavoularis^"*  asymptotic  law 
of  exponential  growth  in  shear  flow.  The  wealth  of  experimental  data  on  homogeneous 
shear  flow  collected  over  the  past  two  decades  appears  to  be  in  agreement  with  this 
hypothetical  physical  picture.  This  is  evidenced  by  Figure  1  containing  several  sets  of 
measured  k  data  versus  non-dimensional  time  that  was  recently  compiled  by  Rohr,  et 
al.i°. 


In  order  for  this  structural  equilibrium  —  with  an  unbounded  exponential  time  growth 
of  k  and  e  —  to  be  valid,  the  higher  order  correlations  in  the  dissipation  rate  transport 
equation  (7)  must  satisfy  certain  consistency  conditions.  When  non-dimensionalized,  (7) 
takes  the  form 


e'  — 


\/Rt- 


2z/  Sxj  dxj 

S 


(17) 


where  Rt  =  A^/i/e  is  the  turbulence  Reynolds  number.  In  deriving  (17),  we  have  made 
use  of  the  fact  that  e  =  i/iviuJi  since  the  turbulence  is  homogeneous.  The  first  of  the  three 
correlations  _  _ 


L1J1OJ2 


dui 
«  J  dxj 


dxj 


dwj 

dxj 


appearing  in  (17)  is  bounded  by  the  Schwarz  inequality.  Direct  numerical  simulations  of 
Rogers^^  for  homogeneous  shear  flow,  as  shown  in  Figure  2,  indicate  that  the  second  of 
these  correlations  (which  is  proportional  to  the  skewness  Sk  in  an  isotropic  turbulence) 
asymptotes  fairly  quickly  to  an  apparent  equilibrium  value  of  0.1.  Hence,  assuming  the 
equilibration  of  ejSk,  it  follows  that  in  order  to  recover  the  exponential  growth  law  (IS) 
for  i*  ^  1,  we  must  have  the  third  of  these  correlations  behave  as 


dwj  duj 
^dxjdx, 

While  (18)  may  appear  to  be  consistent  with  the  traditionally  accepted  scaling  laws  for 
equilibrium  turbulent  flows^®,  it  is  a  somewhat  questionable  possibility  for  homogeneous 
shear  flow  where  one  might  expect  this  correlation  to  equilibrate.  In  particular,  this 
term  represents  the  ratio  of  the  destruction  of  enstrophy  to  the  enstrophy.  The  analo¬ 
gous  quantity  for  the  turbulent  kinetic  energy  —  namely,  the  ratio  of  the  “destruction  of 
turbulent  kinetic  energy”  (e)  to  the  turbulent  kinetic  energy  (A:)  —  does  reach  an  equilib¬ 
rium  value  for  large  times.  Direct  numerical  simulations  of  homogeneous  shear  flow  do 
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not  support  (18),  although  the  issue  is  still  far  from  settled  (see  Figure  3  obtained  from 
Rogers^^,  which  is  not  suggestive  of  a  \/Rt  growth). 

We  will  now  consider  possible  alternative  equilibrium  states  for  homogeneous  shear 
dow.  There  are  two  such  possibilities: 

(1)  an  alternative  structural  equilibrium  where  {e/Sk)oo  —  0,  or 

(2)  a  production-equals-dissipation  equilibrium  with  bounded  energy  states  (i.e.,  with 
koo  <  oo  and  Coo  <  oo). 

The  former  structural  equilibrium  where  (e/Sk)oo  =  0  has  been  shown  by  Speziale 
and  MacGiolla  Mhuiris^^  to  be  primarily  associated  with  solutions  for  k  and  e  that 
undergo  a  power  law  growth  with  time  (i.e.,  for  St  1,  ~  t“,  e  ~  where  a  > 
f3).  However,  these  solutions  are  largely  unstable  within  the  context  of  Reynolds  stress 
transport  models.  Furthermore,  at  the  end  of  all  of  the  previously  conducted  physical  and 
numerical  experiments  on  homogeneous  shear  flow,  e/ Sk  either  appeared  to  equilibrate  to 
a  non-zero  value  or  continued  to  grow  —  results  that  are  not  suggestive  of  an  equilibrium 
state  where  {e/Sk)ao  =  0.  Hence,  we  feel  that  this  alternative  equilibrium  structure  is 
not  a  strong  possibility. 

The  second  possibility  —  namely,  the  production-equals-dissipation  equilibrium  — 
appears  to  contradict  the  physical  and  numerical  experiments  which  indicate  that  the 
turbulent  kinetic  energy  and  dissipation  rate  are  still  growing  at  the  end  of  the  experi¬ 
ments  (i.e.,  for  elapsed  times  St  as  large  as  28).  However,  a  production-equals-dissipation 
equilibrium  wherein 

"^oo  =  Coo  (19) 

(20) 

(with  bounded  values  for  koo,eoo  and  Rt„)  is  consistent  with  the  ensemble  averaged 
Navier-Stokes  equations.  In  the  following  sections,  it  will  be  shown  that  when  the  vortex 
stretching  term  is  accounted  for  in  the  dissipation  rate  transport  equation,  then  this 
production-equals-dissipation  equilibrium  (with  an  early  exponential  time  growth  of  k 
and  e  for  St  <  30)  becomes  a  serious  possibility. 

3.  THE  DISSIPATION  RATE  TRANSPORT  EQUATION  WITH  VORTEX 
STRETCHING 

Batchelor  and  Townsend^®  have  shown  that  the  transient  behavior  of  the  enstrophy 


6 


=  ujiLUi  in  isotropic  turbulence  is  governed  by  the  equation 


where 


'  3  c  3 

3\/r5‘^  ^  3^/TE'^  Rx 


Sk  = 


is  the  skewness  of  the  probability  density  function  of  ~  defined  with  a  negative  sign 

to  make  it  positive  definite.  In  (21),  Rx  =  is  the  turbulence  Reynolds  number 

based  on  the  Taylor  microscale  and  G  is  a  function  defined  by 

G  =  AVr-  (23) 

Here,  Urmj  is  the  root-mean-square  of  the  velocity  fluctuation  component  u,  A  is  the 
Taylor  microscale  derived  from  the  two  point  longitudinal  velocity  correlation  function 
/(r),  and  /’*'  is  the  fourth  derivative  of  /(r)  evaluated  at  r  =  0.  The  first  term  on  the 
right-hand  side  of  (21)  accounts  for  the  effect  of  vortex  stretching  and  is  positive  definite, 
while  the  second  term  is  always  negative  and  leads  to  the  destruction  of  enstrophy. 

Equation  (21)  may  be  converted  into  a  transport  equation  for  e  through  use  of  the 
identity 

e  =  (24) 

which  is  valid  for  any  homogeneous  turbulence.  Since  the  defining  equation  for  the 
microscale  may  be  conveniently  written  as 


it  follows  that  (21)  can  be  converted  to  the  equation 


•  —  ^  I 

3\/r5 


15  k 


Once  values  are  obtained  for  Si^  and  G,  (26)  may  be  solved  together  with  the  kinetic 
energy  equation 

k  =  -t  (27) 

yielding  a  solution  for  isotropic  decay. 
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Batchelor  and  Townsend^®  showed  that  if  G  has  the  form 


(7  = 


(28) 


then  the  solution  of  (26)  -  (27)  is  compatible  with  their  experimental  data  showing  a 
power  law  decay  of  the  kinetic  energy,  with  an  exponent  of  approximately  one,  for  the 
case  of  moderately  large  values  of  Rx.  More  precisely,  the  substitution  of  (28)  into  (26) 
yields  the  equation 


(29) 


which,  when  combined  with  (27),  gives  rise  to  the  exact  solution 


k-ko 


-1 


for  isotropic  decay  (i.e.,  a  power  law  decay  where  k  ~  i“^).  Apart  from  the  specific  value 
of  the  numerical  coefficient  on  the  right-hand  side  of  (29),  which  is  more  commonly  set 
to  a  value  ranging  from  1.83  -  1.92  to  reflect  more  recent  decay  data^®  suggesting  that 
k  ~  or  A:  ~  the  form  of  (29)  is  one  which  forms  a  cornerstone  for  the  standard 
modeled  e  transport  equation  that  is  widely  used  in  turbulence  models. 

In  the  approach  of  Batchelor  and  Townsend,  the  choice  of  G  given  by  (28)  forces  the 
vortex  stretching  term  in  (26)  to  be  exactly  subsumed  by  the  action  of  the  destruction 
of  enstrophy  term.  This  assumption,  as  has  been  noted  above,  guarantees  compatibility 
with  their  isotropic  decay  data.  However,  it  seems  highly  unlikely  that  two  terms  which 
represent  distinct  physical  processes  would  exactly  counterbalance  each  other,  especially 
since  the  destruction  of  enstrophy  term  vanishes  in  the  limit  of  zero  viscosity  whereas 
the  vortex  stretching  term  does  not.  In  fact,  the  survival  of  the  vortex  stretching  term 
in  t!.e  limit  of  zero  viscosity  (i.e.,  the  limit  as  Rt  — >  oo)  is  crucial  for  the  prediction  of 
enstrophy  blow-up  —  a  property  of  solutions  of  the  Euler  equation  (see  Lesieur^°).  This 
can  be  seen  by  setting  u  equal  to  zero  in  (21)  which  yields  the  equation 


dt 
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(0) 

K 


(30) 


where  is  the  zero-viscosity  skewness.  For  constant  >  0,  (30)  predicts  that  the 
enstrophy  blows  up  at  the  critical  time 


tc 


6x/l5  1 
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where  ijj\  is  the  initial  enstrophy.  Although  EDQNM  supports  this  result,  a  finite-time 
enstrophy  blow-up  has  net  been  seen  in  direct  numerical  simulations  of  the  Euler  equation 
(see  Lesieur^°  and  Pumir  and  Siggia^^).  This  means  that  either  is  a  very  small 
constant  or  a  monotonically  decreasing  function  of  time  (in  the  former  case  the  enstrophy 
would  blow-up  at  tc  uJo  whereas  in  the  latter  case  it  would  just  grow  monotonically 
without  bound  as  t  — >  oo).  While  this  issue  has  not  been  fully  resolved,  one  thing  is 
clear:  is  not  identically  zero.  If  =  0,  then  (30)  yields  the  erroneous  prediction 

uj^  =  constant 

which  is  not  supported  by  direct  numerica’  simulations  or  theoretical  analyses  of  the  Euler 
equation  (these  indicate  that  the  enstrophy  grows  dramatically  due  to  vortex  stretching). 

The  standard  modeled  dissipation  -ate  equation  (29)  reduces  to 


in  the  limit  of  zero  viscosity  and,  hence,  incorrectly  predicts  (in  agreement  with  the 
=  0  case)  that  the  enstrophy  is  conserved.  It  is  a  simple  matter  to  show  that  (28)  is 
equivalent  to  the  alternative  scaling 

G  ~  constant,  Sjc  ~  Rx'  (31) 

in  the  e  transport  equation  for  1.  The  scaling  of  G  in  (31)  can  be  obtained  by 

making  the  assumption  that 

dxj  dxj  ^ 

«  - 
(jOkOJk  K 

(namely,  by  the  assumption  that  the  time  scale  of  the  destruction  of  dissipation  is  pro¬ 
portional  to  the  time  scale  of  the  destruction  of  turbulent  kinetic  energy).  While  this 
scaling  for  G  is  quite  reasonable,  the  scaling  of  Sk  in  (31)  is  not  acceptable  since  it  gives 
rise  to  a  vanishing  skewness  in  the  limit  of  infinite  turbulence  Reynolds  numbers  (or  zero 
viscosity).  Hence,  we  will  maintain  the  first  part  of  equation  (31)  but  we  will  not  set 
to  zero.  A  Taylor  series  expansion  of  Sk  in  the  variable  yields 

=  4°^ -f  +  ...  (32) 

where  S'j^^  —  [dSK/dRx^]f^-i_Q.  By  substituting  the  first  part  of  (31)  and  (32)  into  (26), 
we  obtain  the  alternative  modeled  dissipation  rate  equation 

.  ^  4"’  I  „ 

3v/l5  k 
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where  can  be  taken  to  be  a  constant  for  Rx  ^  I  (in  deriving  (33)  we  have  made  use 
of  the  fact  that  Rx^  ~  and  that  j =  R}J^ j k).  In  the  liirut  as  — >  0, 

the  traditional  modeled  transport,  equation  for  e  is  recovered.  We  take  Ce^  =  1.90, 
rather  than  2.0,  since  this  value  yields  a  power  law  decay  in  isotropic  turbulence  where 
k  ~  —  a  value  that  is  more  in  the  range  of  the  most  recent  experiments^®.  As 

alluded  to  earlier,  since  direct  numerical  simulations  of  the  Euler  equation  ha\e  failed  to 
yield  finite  time  enstrophy  blow-up,  the  zero  viscosity  skewness  5^^  is  likely  to  be  small. 
Lagrangian  based  simulations  of  the  Euler  equation  conducted  recently  by  Girimaji  and 
Pope^^  yielded  values  for  5^^  that  were  of  the  order  of  10“^.  Hence,  for  all  of  the 
illustrative  calculations  of  this  paper,  we  will  set 


=  0.01 

for  simplicity  (or  course,  could  more  generally  be  a  monotonically  decreasing  function 
of  time). 

In  figure  4(a),  a  solution  of  (33)  for  isotropic  decay  is  given  for  =  300  (it  should  be 
noted  that  the  new  e-transport  equation  with  vortex  stretching  depends  on  whereas 
the  more  traditional  model  does  not).  It  is  clear  that  for  this  value  of  Ri^,  vortex  stretch¬ 
ing  has  little  effect  on  the  solution.  In  fact,  for  Rt^  <  10,000  —  which  includes  virtually 
all  of  the  physical  and  numerical  experiments  that  have  been  conducted  on  isotropic  tur¬ 
bulence  —  the  differences  are  not  major.  It  is  only  for  extremely  high  turbulence  Reynolds 
numbers  (i.e.,  for  R^  o  >  100,000)  that  vortex  stretching  makes  a  major  difference.  For 
example,  for  Rt„  —  10*",  the  enstrophy  predicted  by  this  new  model  exhibits  a  dramatic 
early  time  growth  which  is  reminiscent  of  the  prelude  to  enstrophy  blow-up  predicted  by 
EDQNM  for  large  values  of  Ri^  (see  figure  4(b)  and  figure  VI  -  4  of  Lesieur^®). 

For  anisotropic  homogeneous  turbulent  flows,  (33)  can  be  generalized  to  the  form 


^  7  SP  ,  ^ 

^  Z^/Ib  ^  k' 


where  is  the  production  of  dissipation  by  mean  strains  given  by  (8).  In  order  to 
achieve  closure,  a  model  for  is  needed.  For  simple  homogeneously  strained  turbulent 
flows,  it  can  be  assumed  that 

oc  V  (35) 


which,  after  invoking  elementary  dimensional  analysis,  yields  the  model 


‘  k  3v/l5  ^  k' 
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where  is  a  dimensionless  constant.  The  model  for  in  (36)  has  been  used  in  the 
k  —  e  model  of  turbulence  as  well  as  in  more  complex  second-order  closures.  Its  success  is 
largely  tied  to  the  fact  that  it  constitutes  a  good  approximation  for  plane  shear  flows^^'^® 
—  the  type  of  flow  being  considered  in  this  study.  For  practical  calculations,  (7^,  can  be 
taken  to  be  1.45  (a  value  obtained  from  equilibrium  shear  flows).  In  the  next  section,  we 
will  apply  this  model  to  homogeneous  shear  flow. 

4.  ILLUSTRATIVE  CALCULATIONS 


In  order  to  illustrate  the  efiect  of  vortex  stretching  on  homogeneous  shear  flow,  we 
will  present  the  results  of  calculations  with  a.  k  -  c  model  for  which  the  Reynolds  stress 
tensor  is  modeled  by 


-C.— 


dv^  dvj 


where  =  0.09  is  a  dimensionless  constant.  While  the  k  —  e  model  is  somewhat 
simplistic  since  it  is  based  on  an  eddy  viscosity,  it  was  recently  shown  by  Speziale  and 
MacGiolla  Mhuiris*®  that  this  model  is  topologically  equivalent  to  the  more  complex 
second-order  closure  models  for  homogeneous  shear  flow  (the  deficiencies  in  the  k  —  e 
model  only  become  pronounced  when  there  are  combinations  of  shear  and  rotation  or 
multi-dimensional  strains).  Hence,  (37)  will  suffice  to  illustrate  the  qualitative  changes 
induced  when  the  effect  of  vortex  stretching  on  the  dissipation  rate  is  accounted  for. 
Equation  (37)  will  be  solved  in  conjunction  with  equations  (6)  and  (36).  The  standard 
k  —  c  model  is  recovered  in  the  limit  as  —>  0. 

For  homogeneous  shear  flow,  the  k  ~  e  model  with  vortex  stretching  yields  the  trans¬ 
port  equations 

k  -  C^—S^  -  e  (38) 

e 


which  are  obtained  by  substituting  (4)  into  (6),  (36)  and  (37). 

For  all  non-zero  values  of  it  is  a  simple  matter  to  show  that  the  solution  to  (38) 
and  (39)  converges  to  an  equilibrium  state,  with  bounded  energies.  The  equilibrium 
values  may  be  found  by  setting  the  right-hand-sides  of  (38)-(39)  to  zero,  yielding  the 
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results 


49  R, 


e^  135  (7^(a,  -  g (SkX 

^o  49  \ej  ■  ^  ^ 

These  relations  clearly  indicate  that  koo/l^o  and  Coo/fo  have  a  Sf^'  dependence  so  that 
the  standard  A:  —  e  model  prediction  of  an  unbounded  growth  of  k  and  e  is  easily  recovered 
in  the  limit  as  5^^  — ^  0. 

Using  (40)  and  (41)  the  following  additional  equilibrium  values  are  also  obtained  for 
this  k  —  e  model  with  vortex  stretching: 


8 

11 

(42) 

(71-^ 

(43) 

135  (a, -aj 

49  q(o)2 

(44) 

These  results  (which  are  independent  of  the  initial  conditions)  differ  from  the  values  of 
(SVe)„  =  (  —uvlk)oo  =  sJo-Cft,  and  Rt^  =  oo  obtained  from  the  standard 

k  —  e  model  where  a  =  (C^j  —  —  1)  «  2.  The  mechanism  by  which  the  presence 

of  the  vortex  stretching  term  has  the  effect  of  creating  bounded  long  time  solutions  lies 
in  its  enhancement  of  the  growth  rate  of  e.  Evidently,  this  increase  in  the  growth  rate 
of  e  is  accompanied  by  a  simultaneous  reduction  in  the  production  of  k,  thus  forcing 
a  production-equals-dissipation  equilibrium.  In  alternative  terms,  vortex  stretching  — 
which  becomes  more  pronounced  at  high  turbulence  Reynolds  numbers  since  it  scales  as 
\/^t  —  eventually  causes  homogeneous  shear  flow  to  undergo  a  saturation  to  an  equilib¬ 
rium  state  with  bounded  component  energies  (the  values  of  which  are  set  by  the  shear 
rate  and  the  viscosity). 

We  will  now  show  that  the  k  ~  t  model  with  vortex  stretching  yields  temporal  evo¬ 
lutions  of  the  turbulence  fields  for  Si  <  30  that  are  in  good  qualitative  agreement  with 
previously  conducted  physical  and  numerical  experiments.  When  non-dimensionalized, 
(38)  -  (39)  take  the  form 


(46) 


e* 


Ce.Cp  — A)*  + 


Co 

SK 


where  again  we  have  =  0.09,  Cej  =  1.45,  Ce,  =  1.90  and  =  0.01  (the  standard 
k  —  e  model  is  obtained  by  setting  =  0).  The  initial  conditions,  which  correspond 
to  an  isotropic  turbulence,  are  taken  to  be  e„/Sko  =  0.296  and  Rt^  =  300.  These  are 
the  approximate  initial  conditions  of  the  large-eddy  simulations  of  Bardina  et  al.^^  which 
will  allow  us  to  make  some  direct  comparisons  between  the  model  and  the  simulations. 

Figures  5(a)  and  6(a)  display  the  short  time  solutions  for  k*  •  d  e*  compared  to  the 
large  eddy  simulation  data.  The  solutions  with  non-zero  vorteA  stretching  are  seen  to 
display  short  term  exponential  growth  in  a  manner  very  similar  to  that  in  the  standard 
k  —  e  closure.  The  effect  of  the  vortex  stretching  term  is  to  reduce  the  growth  rate  of 
k*  and  e*,  though  initially  there  is  a  slight  increase  in  the  magnitude  of  e*.  A  view  of 
these  solutions  over  a  much  longer  time  interval,  as  displayed  in  Figures  5(b)  and  6(b), 
reveals  the  dramatic  effect  that  the  vortex  stretching  term  ultimately  has  on  the  long 
term  growth  of  k*  and  e*.  It  is  seen  that  with  the  vortex  stretching  effect  included,  the 
initial  exponential  growth  rates  are  eventually  suppressed,  so  that  by  St  w  40,  k*  and  e* 
asymf  tote  to  bounded  equilibrium  values. 

Figures  7  and  8  show  the  short  and  long  time  behavior  of  the  dimensionless  ratios  Sk/e 
and  —uvfk  for  the  solutions  obtained  with  and  without  vortex  stretching.  The  curves  in 
Figures  7(a)  and  8(a)  give  the  impression  that  an  equilibrium  state  for  these  quantities 
may  have  been  achieved  by  the  time  St  ~  10.  However,  the  long  time  solutions  reveal 
that,  in  the  case  where  vortex  stretching  effects  are  included,  only  a  local  maximum  is 
reached  —  further  developments  must  occur  before  a  true  equilibrium  is  achieved.  This 
suggests  that  the  apparent  convergence  of  quantities  such  as  Sk/e  seen  in  numerical  and 
experimental  studies,  may  not  signify  that  a  final  equilibrium  state  has  resulted  (i.e.,  it 
may  only  be  a  local  maximum).  The  behavior  of  Sk/e  shown  in  Figure  7(b)  with  vortex 
stretching  present  is  remarkably  consistent  with  the  direct  simulations  of  Rogers  et  al.^^ 
shown  in  Figure  9  (i.e.,  the  tails  in  the  computations  of  Rogers  et  al.^^  may  indicate  that 
an  equilibrium  state  has  not  been  reached). 

Figure  10  provides  a  plot  of  the  long  time  behavior  of  the  computed  turbulence 
Reynolds  number.  It  achieves  an  equilibrium  value  of  approximately  5600,  which  is  more 
than  eighteen  times  its  initial  value.  The  most  significant  effect  of  the  vortex  stretching 
term  on  k*  and  e*,  as  seen  in  Figures  5(b)  and  6(b),  occurs  for  Rt  >  4000.  This  confirms 
the  belief  that  vortex  stretching  is  mostly  a  phenomenon  associated  with  high  turbulence 
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Reynolds  numbers.  Finally,  some  indication  of  the  sensitivity  of  the  computed  solutions 
to  the  skewness  appearing  in  the  vortex  stretching  term  is  shown  in  Figure  11.  This 
contains  the  time  evolution  of  k*  for  a  range  of  values  of  5^^;  as  expected,  increases 
with  decreasing  values  of  Figure  12  shows  the  effect  on  k*  of  a  change  in  the 

initial  values  of  eolSk.  As  would  be  expected  on  physical  grounds,  an  increase  in  the 
dimensionless  shear  rate  leads  to  a  higher  equilibrium  value  for  the  turbulent  kinetic 
energy. 

5.  CONCLUSIONS 

An  alternative  view  concerning  the  equilibrium  structure  of  homogeneous  turbulent 
shear  flow  has  been  presented  based  on  the  inclusion  of  the  effect  of  vortex  stretching.  It 
was  shown  that  the  presence  of  just  a  small  amount  of  vortex  stretching  can  ultimately 
lead  to  a  production-equals-dissipation  equilibrium  with  bounded  energy  states.  For 
elapsed  times  <  30  —  which  includes  the  largest  values  of  St  considered  in  any  of  the 
previously  conducted  physical  and  numerical  experiments  —  the  introduction  of  a  small 
vortex  stretching  term  into  the  standard  modeled  dissipation  rate  transport  equation  still 
yields  an  exponentially  growing  kinetic  energy  and  dissipation  for  early  times.  However, 
since  this  vortex  stretching  term  scales  as  %/7^,  it  eventually  becomes  dominant,  causing 
a  saturation  of  the  system  to  an  equilibrium  state  with  bounded  turbulent  kinetic  energy 
and  dissipation.  Although  this  alternative  physical  picture  of  homogeneous  shear  flow  is 
contrary  to  the  commonly  accepted  asymptotic  laws,  for  which  an  unbounded  exponential 
time  growth  of  k  and  e  is  postulated,  it  is  a  real  possibility  that  needs  to  be  seriously 
considered  in  the  future.  New  physical  and  numerical  experiments,  for  larger  elapsed 
times  St,  could  shed  more  light  on  the  issue.  The  recent  experiments  of  Tavoularis 
and  Karnik^,  which  were  conducted  up  to  St  =  28,  did  show  some  tendency  of  the 
integral  length  scales  to  level  off  —  a  feature  which,  if  more  solidly  established,  would 
be  supportive  of  the  existence  of  a  production-equals-dissipation  equilibrium.  However, 
it  is  possible  for  the  kinetic  energy  and  dissipation  rate  to  equilibrate  while  the  integral 
length  scales  grow  without  bound  (perhaps  this  is  the  more  likely  possibility  since  the 
flow  field  is  infinite). 

Finally,  some  comments  should  be  made  concerning  the  implications  that  the  results 
of  this  paper  have  on  turbulence  modeling.  Since  all  of  the  commonly  used  two-equation 
models  and  second-order  closures  based  on  the  turbulent  dissipation  rate  neglect  vortex 
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stretching,  they  predict  an  unbounded  exponential  time  growth  of  k  and  e  in  homo¬ 
geneous  shear  flow.  This  type  of  behavior  has  been  shown  to  cause  problems  in  the 
calculation  of  certain  inhomogeneous  turbulent  flows.  The  singularity  in  plane  stagna¬ 
tion  point  turbulent  flow  represents  a  prime  example^^.  Hence,  if  it  is  ultimately  justified, 
the  alteration  of  turbulence  models  to  yield  a  production-equals-dissipation  equihbrium 
in  homogeneous  shear  flow  via  vortex  stretching  could  make  their  behavior  more  robust 
in  other  turbulent  flows.  A  full  resolution  of  this  issue  is  only  likely  to  come  from  a  rig¬ 
orous  mathematical  proof  based  on  some  appropriate  energy  method.  For  the  meantime, 
however,  the  results  of  this  study  clearly  establish  the  need  to  re-examine  this  issue. 
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Figure  1.  Experimental  data  on  homogeneous  shear  flow  measured  b; 
Corrsin®  and  Tavoularis  and  Karnik®  (as  compiled  by  Rohr  et  al.*°). 


Figure  4.  Comparison  of  solutions  for  isotropic  decay:  - =  0.01;  -  -  -  -  standard 

k  —  e  model,  (a)  Time  evolution  of  the  turbulent  kinetic  energy  for  Rt^  =  300,  and  (b) 
time  evolution  of  the  enstrophy  for  iZ*.  =  10®. 


Figure  4.  Comparison  of  solutions  for  isotropic  decay;  - =  0.01;  -  -  -  -  standard 

k  —  €  model,  (a)  Time  evolution  of  the  turbulent  kinetic  energy  for  Rt^  =  300,  amd  (b) 
time  evolution  of  the  enstrophy  for  Rt.  =  10®. 


Figure  5.  Comparison  of  the  model  predictions  with  the  large-eddy  simulations  of  Bar¬ 
dina,  et  al.^^:  o  k*  from  large  eddy  simulation^^; - k*  obtained  with  =  0.01;  -  -  - 

-  k*  obtained  from  the  standard  k  —  e  model,  (a)  Short  time  solution,  and  (b)  long  time 
solution. 


Figure  5.  Comparison  of  the  model  predictions  with  the  large-eddy  simulations  of  Bar¬ 
dina,  et  al.^*;  o  k*  from  large  eddy  simulation'^; - k*  obtained  with  =  0.01;  -  -  - 

-  k*  obtained  from  the  standard  k  —  €  model,  (a)  Short  time  solution,  and  (b)  long  time 
solution. 
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Figure  6.  Comparison  of  the  model  predictions  with  the  large-eddy  simulations  of  Bar¬ 
dina,  et  al.^^;  o  e*  from  large  eddy  simulation^^; - e*  obtained  with  SP  =  0.01; - 

-  e  obtained  from  the  standard  k  —  e  model,  (a)  Short  time  solution,  and  (b)  long  time 
solution. 


0 


20 


40 

St 


60 


80 


Figure  6.  Comparison  of  the  model  predictions  with  the  large-eddy  simulations  of  Bar¬ 
dina,  et  al.'*:  o  e*  from  large  eddy  simulation**; - e*  obtained  with  =  0.01;  -  -  - 

-  e*  obtained  from  the  standard  k  —  e  model,  (a)  Short  time  solution,  and  (b)  long  time 
solution. 
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Figure  11.  Sensitivity  of  the  model  predictions  for  the  turbulent  kinetic  energy  to  (a) 
SP  =  0.0133;  (b)  SP  =  0.01;  (c)  sP  =  0.0066;  (d)  sP  =  0.005  and  (e)  sP  =  0.0033. 
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